Submitted to ApJ 2007 May 12 

Preprint typeset using I^TgX style emulateapj v. 1 1/26/04 



X-RAY EMISSION FROM PSR J1809-1917 AND ITS PULSAR WIND NEBULA, POSSIBLY ASSOCIATED WITH THE 

TEV GAMMA-RAY SOURCE HESS J1809-193 

O. Kargaltsev and G. G. Pavlov 

The Pennsylvania State University, 525 Davey Lab, University Parle, PA 16802, USA 
Submitted to ApJ 2007 May 12 



o 
o 



o 



> 
oo 
l> 
m 

in 
o 
i> 
o 



X 



ABSTRACT 

We detected X-ray emission from the 50-kyr-old pulsar J1809-1917 and resolved its pulsar wind nebula 
(PWN) with the Chandra X-ray Observatory. The pulsar's observed flux is Fps^ = (1.8 ±0.2) x 10"'"* ergs 
cm"^ s"' in the 1-6 keV band. A two-component blackbody-npower-law (BBh-PL) fit of the pulsar's spectrum 
yields the photon index Fpsr = 1.2 ±0.6 and luminosity Lpsr = (4± 1) x 10^' ergs s"' of the PL component, 
in the 0.5-8 keV band, for a plausible distance d = 3.5 kpc and nu = 0.7 x 10^^ cm~^. The BB component 
corresponds to the temperature T « 2 MK, and bolometric luminosity Lboi 1 x lO-'^ ergs s"'. The bright 
inner PWN component of a 3" x 12" size is elongated in the north-south direction, with the pulsar close to its 
south end. This component is immersed in a larger (« 20" x 40"), similarly elongated outer PWN component 
of lower surface brightness. The elongated shape of the compact PWN can be explained by the ram pressure 
confinement of the pulsar wind due to the supersonic motion of the pulsar The observed flux of the compact 
PWN, including both components, is Fp^„ ~ (L5 ±0.1) x 10"'^ ergs cm"^ s~' in the 1-6 keV band. The 
PWN spectrum can be fitted with a PL model with tin « 0.7 x 10^^ cm"^ and photon index Fpwn = 1.4±0.1, 
corresponding to the 0.5-8 keV luminosity Lp^n « 4 x 10^^ ergs s"'. The compact PWN appears to be inside 
a more extended (« 4' x 4') emission with the total observed flux Fgxt ~ 5 x 10"'"' ergs s~' in the 0.8-7 keV 
band. This large-scale emission is more extended to the south of the pulsar, i.e. in the direction of the alleged 
pulsar motion. To explain the extended X-ray emission ahead of the moving pulsar, one has to invoke strong 
intrinsic anisotropy of the pulsar wind or assume that this emission comes from a relic PWN swept by the 
asymmetrical reverse SNR shock. The pulsar and its PWN are located within the extent of the unidentified TeV 
source HESS J1809-193. The brightest part of the TeV source is offset by ~ 8' to the south of the pulsar, i.e. in 
the same direction as the large-scale X-ray emission. Although the association between the PSR J1809-1917 
and HESS J1809-193 is plausible, an alternative source of relativistic electrons powering HESS J1809-193 
might be the serendipitously discovered X-ray source CXOU J 180940.7- 192544. In addition to the CMBR 
or Galactic starlight background, the low-frequency seed photons for Compton upscattering to TeV energies 
might be supplied by bright infrared emission from dust-molecular clouds seen within HESS J1809-193. 
Subject headings: pulsai's: individual (PSR J1809-1917) — X-rays: individual (CXOU J180940.7-192544, 
CXOU 180933.3-192959) — gamma-rays: individual (HESS J1809-193) — ISM: indi- 
vidual (IRAS 18067-1927, IRAS 18067-1921) 



1. INTRODUCTION 

Chandra and XMM-Newton observations have established 
the ubiquity of X-ray pulsar wind nebulae (PWNe) around 
young rotation-powered pulsars (see the reviews by Kaspi et 
al. 2006 and Gaensler & Slane 2006). The X-ray PWN emis- 
sion is produced by relativistic particles gyrating in the mag- 
netic field downstream of the termination shock in the pulsar 
wind (Kennel & Coroniti 1994; Arons 2004). Most of the 
PWNe have been discovered around young (t < 30 kyr) pul- 
sars. The innermost parts of the young PWNe often show 
axisymmetric morphologies, including toroidal structures and 
jets along the pulsar's spin axis. Recently, it has become ap- 
parent that PWNe accompanying older pulsars can also be 
quite luminous (e.g., McGowan et al. 2006). Many of these 
older PWNe exhibit cometary morphologies indicating that 
the pulsar wind is confined by the ram pressure caused by 
the supersonic motion of the pulsar in the ambient medium. 
Studying X-ray PWNe of various ages helps understand the 
nature and evolution of the ultrarelativistic pulsar winds and 
their interaction with the ambient medium. 

An interesting object for such investigations is PSR J809- 
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1917 (hereafter J1809). The discovery of this radio pulsar 
(P = 82.7 ms) in the Parkes Multibeam Pulsar Survey' was 
reported by Morris et al. (2002). The pulsar's dispersion mea- 
sure, DM = 197cm~^pc, and the Galactic electron distribu- 
tion models by Taylor & Cordes (1993) and Cordes & Lazio 
(2002) give the distance to the pulsar of 3.7 and 3.5 kpc, re- 
spectively. Having the spin-down age r=P/2P = 51 kyr and 
spin-down power £ EE 47r2/Pf -3 ~ 1.8 x lO^*" ergs s"', J 1809 
is somewhat older and less energetic than the famous Vela pul- 
sar (r = 1 1 kyr, £ = 6.9 x lO-'^ ergs s"'), which is accompanied 
by a remarkable PWN resolved in radio (Dodson et al. 2003b) 
and X-rays (Pavlov et al. 2003, and references therein). How- 
ever, it is much more energetic than typical "middle-aged" 
pulsars, such as B0656H-14 (r = 1 10 kyr, £ = 3.8 x 10^"* ergs 
s"') and Geminga (r = 340 kyr, £ = 3.2 x 10^"* ergs s"'), whose 
PWNe are very faint (e.g., Pavlov et al. 2006). 

J 1809 is young enough to look for a remnant of the super- 
nova that created the pulsar Deep radio observations by Bro- 
gan et al. (2004) have revealed two compact HII regions and 
two SNRs, G 11.03-0.05 and G11.18h-0.11, projected near 
the pulsar (the pulsar's Galactic coordinates are / = 11.094°, 

' http://www.atnf.csiro.au/reseai'ch/pulsar/pmsurv 
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b = +0.080°). The distances to the SNRs, estimated from 
the radio surface brightness-diameter (S-Z)) relation (Case 
& Bhattacharya 1998), are ~ 16 and ~ 17 kpc, respectively. 
However, given the very large uncertainty of the S-Z) rela- 
tion for faint SNRs, the association of J1809 with one of the 
SNRs cannot be ruled out despite the discrepant distance es- 
timates. Brogan et al. (2004) estimated that the pulsar would 
need to have the transverse speed of about 200 or 140 km s~' 
to originate from the geometrical center of G 11.03-0.05 or 
Gl 1.18+0. 11, respectively, assuming that all three are at the 
same distance of 4 kpc and the pulsar's true age is 50 kyrs. 
The estimated speeds are close to the average speeds recently 
measured for a large sample of radio pulsars by Hobbs et al. 
(2005). 

Before the detailed radio studies were carried out, the re- 
gion had been observed by ASCA, first as part of the Galac- 
tic plane survey (Sugizaki et al. 2001), and then with a 
deeper follow-up exposure (Bamba et al. 2003; Ueno et al. 
2005). These observations revealed an amorphous diffuse 
emission that encompassed the pulsar position. No pulsar 
has been detected in these X-ray observations. The observed 
large-scale emission was attributed to a new SNR, dubbed 
Gl 1.0+0.0. Bamba et al. (2003) have found that the X-ray 
spectrum of the putative SNR fits the absorbed power-law 
(PL) model with photon index F = 1.6;!^ hydrogen column 
density nH,22 = "h/CIO^^ cm'^) = 0.8 ± 0.3, and 0.7-10 keV 
flux of 3.8 X 10"'^ ergs cm"^ s"'. Using the best-fit «h and 
assuming that the mean density in the Galactic plane is 1 H 
cm~^, Bamba et al. (2003) estimated the distance of 2.6 kpc, 
which gives the X-ray luminosity of ss 3.7 x lO^"' ergs s"' (in 
0.7-10 keV). These authors suggested that Gl 1.0+0.0 could 
be a Crab-like (plerionic) SNR. 

Studying of the J 1809 field has become particularly in- 
teresting after the recent discovery of the TeV 7-ray source 
HESS J1809-193 (hereafter HESSJ1809; Aharonian et al. 
2007). The brightest, firmly detected part of this extended 
source (radius 12-15') is centered at R.A.= 18*'09.8'", 
decl.= -19°25', about 8' from J1809; its 7-ray flux is w 
1.4 X 10"^^ ergs cm"^ s~^ in the 1-10 TeV band. Given the 
small angular separation, it seems plausible that the HESS 
source could be powered by this pulsar, in which case its 7- 
ray luminosity, L^k,2x lO^dl^ ergs s"\ would be about 1% 
of the pulsar's spindown power. 

In this paper, we describe the results of a Chandra obser- 
vation of PSR J1809-1917 and its compact synchrotron neb- 
ula^ and discuss its possible connection to HESS J 1809. We 
also describe the multiwavelength properties of objects lo- 
cated within the central part of HESS J 1809, including two 
newly discovered X-ray sources, and discuss their relation to 
HESS J1809. The details of the observation and the data anal- 
ysis are presented in §2. In §3 we discuss possible interpre- 
tations of the PWN morphology, describe inferences from the 
pulsar spectrum, and speculate on the nature of HESS J 1809 
and its relation to the other sources in the fleld. Our main 
results are summarized in §4. 

2. OBSERVATIONS AND DATA ANALYSIS 

J 1809 was observed with the Advanced CCD Imaging 
Spectrometer (ACIS) on board Chandra on 2004 July 21 (Ob- 

sID 3853). The useful scientific exposure time was 19,955 s. 
The observation was carried out in Faint mode, and the pul- 

^ Preliminary results of this observation have been presented by Sanwal et 
al. (2005). 



sar was imaged on S3 chip, « 0.72' off-axis. The other ACIS 
chips activated during this observation were SI, S2, S4, 12, 
and 13. The detector was operated in Full Frame mode, which 
provides time resolution of 3.24 seconds. The data were re- 
duced using the Chandra Interactive Analysis of Observations 
(CIAO) software (ver. 3.2.1; CALDB ver. 3.0.3). 

2.1. Images 

Figure 1 shows the ACIS-S3 image of the region around 
J1809. An extended X-ray source is clearly seen in the image 
around R.A.= 18'>09'"43.123'\ decl.= -19°17'38.17" (these 
are the coordinates of the center of the brightest pixel). The 
difference of 0.5" between this position and the radio position 
from Morris et al. (2006) is within the uncertainty of absolute 
Chandra astrometry (0.6" at the 90% confidence level). The 
close match between the X-ray and radio positions and the ex- 
tended morphology of the observed X-ray emission allow us 
to conclude that we detected the X-ray emission from J 1809 
and its PWN. 

The brighter inner PWN component of an w 12" x 3" size 
(i.e. 0.2 X 0.05 pc^ at <i = 3.5 kpc) is elongated along the ap- 
proximate symmetry axis (position angle « 14° east of north). 
The linear profile of the surface brightness distribution along 
the symmetry axis, extracted with the 4.9" x 0.49" (i.e. 10x1 
pixels) sliding box, is shown in Figure 2. The inner PWN is 
surrounded by a similarly elongated ~ 20" x 40" "halo" of 
lower surface brightness, which we will call the outer PWN 
component (see Fig. 1, top). In addition to these two relatively 
compact components, the heavily binned and smoothed im- 
age (see Fig. 1, bottom right, and Fig. 3a) reveals even fainter 
large-scale extended emission, concentrated at the lower half 
of the S3 chip (south of the pulsar). The morphology of this 
emission (see also the ASCA image in Fig. 3c) indicates that 
it possibly extends further south but may not be discernible on 
the S2 chip because of its lower sensitivity. 

The only source on the S2 chip detected above the 3cr level 
is located at R.A. = 18>'09'"40.725^ decl. = -19°25'44.10" 
(the Icr centroid uncertainty is 0.28" and 0.29" in right as- 
cension and declination, respectively), well within the bright- 
est central part of HESS J 1809 (about 2' west of the peak of 
the TeV brightness distribution; see Fig. 3). We designate 
this source CXOU J180940.7-192544 and call it Chi here- 
after Although Chi appears to be extended in the ACIS im- 
age, the point spread function (PSF) simulation shows that 
this is Ukely the result of the off-axis location (off-axis angle 
6 = I'A). We find no significant large-scale non-uniformities 
in the X-ray background on the S2 chip (see Fig. 3). 

We have also examined the SI chip image and found 
only one source detected above the 3(t level. The source, 
CXOU 180933.3-192959 (hereafter Ch2), is located at 
R.A. = 18'>09■"33.336^ decl. = -19°29'59.89" (tiie Itr cen- 
troid uncertainty is 0.70" in R.A. and 0.74" in decl.), which 
is about 6' from the center of HESS J 1809. The X-ray 
source is consistent with being point-like; however, due to 
the broadened PSF (FWHM « 10") at flie large off-axis an- 
gle (11.8'), a compact extended source cannot be ruled out. 
We have also searched for diffuse emission features on the 
SI chip and produced an image corrected for the exposure 
map non-uniformities (see Fig. T>b). In this image the X-ray 
emission is systematically brighter toward the western edge 
of the SI chip, likely due to the imperfection of the mono- 
energetic (2 keV) exposure map correction at large off-axis 
angles (see, however, §2.5). We found no traces of the ra- 
dio SNRs Gl 1 .03-0.05 and Gl 1 . 1 8+0. 1 1 in the ACIS images 
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Fig. 1.— Top left: ACIS-S3 image of J1809 and its PWN (0.8-7 keV; pixel size 0.49"). Top right: Extraction regions used for tlie analysis of the PWN 
components (0.8-7 keV; pixel size 0.98"). Bottom left: Adaptively smoothed sub-pixel resolution image (0.8-7 keV; pixel size 0.25") obtained by removing 
the pipeline pixel-randomization and applying the sub-pixel resolution tool (based on analyzing the charge distribution produced by an X-ray event; Tsunemi et 
al. 2001; Mori et al. 2001). The straight line shows the approximate symmetry axis of the X-ray PWN. Bottom right: Heavily binned (pixels size 3.94") and 
smoothed (with a gaussian kernel of r = 25") ACIS-S3 image of J1809 and its PWN. The brightness and smoothing scales are chosen to show the fainter, more 
extended emission. 



(see Fig. 3). 



2.2. Spectral analysis 
2.2.1. PWN spectrum 



spectra for each of the two PWN regions with the absorbed 
PL model, allowing the hydrogen column density, nu, to 



vary. These fits result in spectral slopes Fp 



1.2-2.0 and 



We extracted the PWN spectra from two elliptical regions 
shown in Figure 1 (top right). The smaller elliptical region 
(of 25.1 arcsec^ area) encompasses the brighter inner PWN, 
while the larger elliptical region (of 399.4 arcsec^ area) in- 
cludes the outer PWN component of a lower surface bright- 
ness. To avoid the contamination of the PWN spectrum by 
the pulsar, we excluded from these regions the circular region 
of 1.46" radius centered on the brightest pixel. The back- 
ground (367 counts in 5202 arcsec" area, 0.3-8 keV band) 
was measured from the 37" < r < 55" annulus centered on 
the pulsar The total numbers of counts extracted from the 
smaller and larger elliptical regions (excluding the 1 .46" ra- 
dius circle) are 153 and 365, of which 99.2% and 92.4 % are 
expected to come from the source, which gives 151 ± 12 and 
337 ± 19 PWN counts in the two regions. The observed PWN 
fluxes (in the 1-6 keV band) are Fi„ner = (8.1 ± 0.7) x lO"'"* 
and Fpwn = (14.7 ±0.8) x 10"'"^ ergs s"' cm"^ for the inner 
and the entire (inner-i-outer) PWN, respectively. The corre- 
sponding average intensities are /inner = (4.4 ±0.4) x 10"'^ and 
/p„n = (3.7 ±0.2) X 10"'^ ergs cm"^ s"' arcsec"l 

To investigate the PWN spectral properties, we first fit the 



«H,22 ~ 0.5-1 . 1 (the ranges correspond to the 68% confidence 
level for a single interesting parameter). The difference in 
the best-fit parameters for the inner and the entire PWN (see 
Fig. 4) is statistically insignificant. In particular, we see no 
spectral softening (expected due to synchrotron cooling) in 
the spectrum extracted from the larger region (entire PWN) 
compared to the spectrum of the inner PWN^. Therefore, be- 
low we will use the better constrained best-fit parameters for 
the entire PWN (see Figs. 5 and 6) whenever we refer to the 
PWN spectral properties. 

Given the J1809's dispersion measure, DM = 197 cm"-^ 
pc (i.e., electron column density = 6.08 x 10^*' cm"^), the 
nu,22 value of 0.72 (obtained from the fit to the spectrum of 
the entire PWN) corresponds to the ISM ionization degree 
Hg/nu ~ 8.4%, only slightly below the commonly used value 
of 10%, which corresponds to nu.22 = 0.61. The total Galactic 
HI column density in this direction is 1 .8 x 10^^ cm"^ (Dickey 

^ We have also measured the spectrum of the outer PWN separately and 
found that the best-fit PL parameters are consistent with those obtained for 
the inner PWN and the entire PWN spectra, but their uncertainties are larger 
because of the larger background contribution. 
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Fig. 2. — Top: One-dimensional surface brightness distribution along the 
symmetry axis of the J1809 PWN. Bottom: Extraction regions used to mea- 
sure the pulsar's spectrum (see §2.2.2) and the sliding box used for measuring 
the one-dimensional surface brightness distribution (see §2.1 for details). 

& Lockman, 1990). 

The isotropic luminosity of the entire (inner-i-outer) X-ray 
PWN is Lp„„ EE 47r£/2Fp™f^ « (3.9 ±0.3) x lO^^^^^ ergs s"' 
in the 0.5-8 keV band (1^3 5 = (i/3.5kpc), with approximately 
equal contributions from the outer and the inner PWN com- 
ponents (see Table 1). 

We have also attempted to fit the spectrum of the faint 
large-scale emission (within the ellipse shown in Fig. 3a, 
excluding the compact PWN; 13.2 arcmin^ area) surround- 
ing J 1809 and its compact PWN. Because of the low sur- 
face brightness of this emission and large background con- 
tribution (~ 65%), the spectral fits yield inconclusive results. 
With the absorption being fixed at nH,22 = 0.72, the fits with 
single PL and thermal plasma"* models give F = 2.3 ± 0.3 
(xl « 2) and kT = 0.67 ± 0.08 keV, (xl « 3). The main con- 
tribution to the large values of xi comes from high-energy 
channels (> 4 keV and > 2 keV for PL and mekal, respec- 
tively), which suggests either a mixture of thermal and hard 
non-thermal emission or varying F or kT within the extrac- 
tion region. The measured flux, independent of the model, 
is f'ext '~ 5 X 10"'"' ergs s~' cm~^ in the 0.8-7 keV band, cor- 
responding to the average surface brightness 4x1 1 x 10"'^ 
ergs cm"- s"' arcsec"^. The unabsorbed flux, obtained from 
the PL model, is F™"*"* ~ 8 x 10"'^ ergs s"' cm"^ in the 0.5-8 
keV band. 

^ The model "mekal" in XSPEC, with standard abundance 



TABLE 1 
PL FITS TO THE PWN SPECTRUM 



Model 


"H,22 


AT" 


r 


(C or x^)''/dof 




Entire PWN 
Inner PWN 


0.72 
0.98 




1 41+0.12 


0.82/34 
335/526 


-, 71+0.19 
^■^^-0.17 



Note. — The fits are for fixed «h.22 = "h/IO"^ cm ^ . The uncertainties 
are given at 68% confidence level for a single interesting parameter. 

"Spectral flux in units of 10"^ photons cm"- s"' keV"' at 1 keV. 

''We use the C statistic (Cash 1979) for the inner PWN (which has less 
counts) and the statistic for the entire PWN. 

'^Unabsorbed isotropic luminosity in the 0.5-8 keV band, in units of 10^^^ 
ergs s"' . 

2.2.2. Pulsar spectrum 

To minimize the contamination by the PWN, the pulsar 
spectrum was extracted from a small circular aperture (green 
circle in the bottom panel of Fig. 2) with the radius of 1.5 
ACIS pixels (~ 0.74", 85% encircled energy radius), while 
the background was taken from the 10 arcsec^ region between 
the white circle and white ellipse in the bottom panel of Fig. 2. 
The background region includes the bright part of the PWN; 
it contributes w 8 counts to the total of 67 counts extracted 
from the source aperture. Given the small number of counts 
and the large background contribution, we chose not to sub- 
tract the background but rather to fit it simultaneously with 
the source spectrum, using an additional absorbed PL model 
with the same hh as for the source. The pulsar's absorbed flux 
isFp.sr = (1.8±0.2)x 10"''*ergs cm"^ s"' in the l-6keVband 
(aperture corrected). 

Although the absorbed PL model formally fits the pulsar 
spectrum indicating a soft PL (Fps, « 2.6; see Table 2 and 
Figs. 7 and 8), the fit yields nH.22 ~ 0.4, smaller than that 
for the PL fit to the entire PWN spectrum. To obtain a bet- 
ter constrained fit, we fixed the hydrogen column density at 
nH.22 = 0.72, obtained above from the PL fit to the PWN spec- 
trum. With this nil, the single PL fit is still acceptable, but 
it yields a large photon index, F = 3.2 ±0.4, suggesting a 
thermal emission contribution. A two-component, BB-i-PL, 
fit yields reasonable values of fitting parameters (see Table 2 
and Figs. 7-9), which, however, are poorly constrained be- 
cause of the small number of photons detected. The slope of 
the PL component is Fpsr = 1 .2 ± 0.6, and its unabsorbed lu- 
minosity is Lpsi- ~ 4 x 10""(i| 5 ergs s"' in the 0.5-8 keV band. 
The temperature and the projected area of the BB component 
are strongly correlated (see Fig. 9), which results in large un- 
certainties for these parameters. The best-fit temperature for 
the BB component is T « 2 MK, while the projected emit- 
ting area, ^ ^ 3 x lO^iifj m^, is smaller than that of the sur- 
face of a neutron star (tt/?^ 3 x 10^ m^), but larger than the 
conventional polar cap area Ape = lii^R^ /cP sa 2.5 x 10^ m^. 
The corresponding bolometric luminosity, Lboi = ^ 
1 X 10^2^15 ergs s"i. 

2.2.3. Chi spectrum 

We extracted 99 events from the elliptical region (with the 
semi-major and semi-minor axes of 10.3" and 6.3", respec- 
tively) centered on the best-fit position of Chi and conform- 
ing to the shape of the off-axis PSF. Based on the PSF simu- 
lation, the extraction region contains « 95% of point source 
counts. The expected background contribution, 4.4 counts, is 
negligible (w 4% of the total number of counts). The mea- 
sured source flux in 0.8-7 keV is /thi = (8.6 ±0.9) x 10"'"* 
ergs cm"^ s"' , after correcting for vignetting and finite extrac- 
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Fig. 3. — J1809 and its vicinity at different wavelengtlis. All the eight panels show the same area on the sky. a: ACIS-S3, S2 and SI images (0.8—7 keV; 
smoothed with 23.6" gaussian kernel). The white ellipse shows the region used to estimate the flux of the large-scale diffuse emission (see §2.2.1). h: The same 
image divided by the exposure map with point sources removed (except for J1809). The contours show the radio emission from two SNRs and two compact 
(possibly HII) regions (labeled A and B; adopted from Brogan et al. 2004). c: ASCA GIS2 and GIS3 combined image (39 ks total exposure; 0.5—10 keV) of 
the same region (the image has been divided by the exposure map; astrometry has been corrected, courtesy of E. Gotthelf). d: Combined 9 ks EPIC-PN and 
MOSl+2 image (0.5-10 keV) obtained in the Galactic plane survey (PI: R. Warwick), e: The same as in the panel a with no smoothing appUed (pixel size is 
7.4"). /■ NRAO VLA Sky Survey (NVSS; Condon et al. 1998) image at 1.4 GHz. g: Spitzer IRAC 8 ^im image from the GLIMPSE survey, h: Spitzer MIPS 
24 fim image. The white circle in panels c. f, g, and h marks the position of the pulsar, the diamond and the star mark the Chi and Ch2 positions, respectively, 
and the cross shows the position of the peak of the TeV brightness distribution of HESS J1809. 

tion aperture size. The spectrum fits well by an absorbed PL 
model (see Table 3) with «h.22 = 1 .2 ± 0.4 and T = 1 .4 ± 0.4. 
The unabsorbed flux in 0.5-8 keV is F™f ^ « 1 .4 x lO"'^ ergs 

-2 -I 

cm s . 

2.2.4. Ch2 spectrum 

We extracted 90 events within a 30" radius around the best- 
fit position of Ch2. Based on the PSF simulation, the extrac- 
tion region contains 85% of the source counts if Ch2 is a 
point source. The background contributes « 32 counts, (i.e. 
« 36% of the total number of counts). The measured 0.8-7 
keV source flux is Fchi ~ 9 x 10"'"^ ergs cm"^ s"', after cor- 
recting for vignetting and finite extraction aperture size (as- 
suming a point source). The thermal plasma (mekal) model 
fits the Ch2 spectrum best, although fits with absorbed PL and 
BB models are also acceptable (see Table 3 for details). 



TABLE 2 
Fits to the pulsar spectrum 



Model »H.22 AA" or A'' T or kT^ C Lx** or Lboi' 



PL 


0.72 


22.6^-^ 


T 99-K).42 
•'■^ -0.38 


412 


0.98-:- 


PL 


0.4 




9 sQ+fl.34 
^■-^ -0.28 


411 


0.65!«;>« 


PL+BB(PL) 


0.72 


2 04+2.156 


1.23 ±0.62 


407 


3T°-^^ 

"■-^'-0.10 


PL+BB(BB) 


0.72 




"■^ -0.03 


407 


1 tt*^-^ 

^•"-0.4 



Note. — The fits are for fixed nH,22 = "h/IO"^ cm"^. The uncertainties 
are given at 68% confidence level for a single interesting parameter. In each 
case the fits were done using the C statistics (Cash 1979) and the unbinned 
source and background spectra, with the total of 1052 channels. 

"Spectral flux in units of 10"* photons cm"- s keV"' at I keV. 

''Projected area of the emitting region for the BB model, in units of lO*" m-. 

■^BB temperature in keV. 

''Unabsorbed PL luminosity in the 0.5-8 keV band, in units of 10^^ ergs s"' . 
■^Bolometric BB luminosity, in units of 10^- ergs s"'. 
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Fig. 4. — Confidence contours (68% and 90%) in the nn-F plane for the 
PL fit to the inner (dashed) and entire (solid) PWN spectra. The contours are 
obtained with the PL normahzation fitted at each point of the grid. 
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Fig. 5. — Entire (outer+inner) PWN spectrum fitted with the PL model (see 
Table 1 for details). 
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Fig. 6.— Confidence contours (68% and 90%) for the PL fit to the entire 
PWN spectrum (for a fixed nH,22 = 0.72). The PL normalization is in units 
of 10~* photons cm"^ s~' keV~' at 1 keV. The dashed curves are the loci 
of constant unabsorbed flux in the 0.5-8 keV band; the flux values near the 
curves are in units of lO"''* ergs cm~^ s~^ 




Fig. 7. — Top: Pulsar spectrum fitted with the absorbed PL model. Middle: 
Pulsar spectrum fitted with the absorbed PL model where the nn was fixed 
at the best-fit value for the entire PWN (nH,22 = 0.72). Bottom: The pulsar 
spectrum fitted with the PL+BB model with fixed nH,22 = 72. The dashed 
and dash-dotted histograms correspond to the BB and PL components, re- 
spectively. The residual panels show the contributions of the energy bins into 
the best-fit C-statistic (multiphed by -1 when the niraiber of data counts is 
smaller than the number of model counts). 

2.3. Timing 

The 3.24 s time-resolution of this observation is insufficient 
to observe the 82.7 ms pulsations from the pulsar. We found 
no significant variability in the Chi and Ch2 Ughtcurves. 

2.4. Archival X-ray data 

The region of interest has been previously observed by 
ASCA for 39 ks (Bamba et al. 2003). The ASCA GIS im- 
age in Figure 3c shows a region of enhanced X-ray bright- 
ness that encompasses J 1809 and its compact PWN; how- 
ever, the diffuse X-ray emission also extends at least ~ 10' 
southward, covering the central region of the brightest piirt of 
HESS J1809. 

The field was also partly observed by the XMM-Newton 
as a part of the Galactic plane survey (PI: R. Warwick). A 
point-like object is clearly seen at the pulsar position in the 
combined EPIC (MOS1-hMOS2-hPN) image shown in Figure 
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Fig. 8. — Confidence contours (68% and 90%) for the PL fit to the pulsars's 
spectrum with «h.22 = 0.72 (green), PL+BB fit with nH,22 = 0.72 (red) and PL 
fit with JiH,22 = 0.4 (black) which is the best-fit value for the PL model. The 
PL normalization is in units of 10"^ photons cm"^ s"' keV"' at 1 keV. The 
dashed curves are the lines of constant unabsorbed flux in the 0.5—8 keV band 
(the flux values are in units of lO"''' ergs cm"^ s"'). 
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BB+PL fit to the pulsar's spectrum, for fin, 22 = 0.72. The BB normaHzation 
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distance of 3.5 kpc. The fines of constant bolometric flux (in units of lO"'** 
ergs cm"^ s"'), are plotted as dashed lines. 



TABLE 3 
Fits to the ChI and Ch2 spectra 



Source Model 


«H.22 


AT" or Ji^ 


r or kT" 


(Corx')''/dofF7orFbol'= 


Chi PL 
Ch2 PL 


i.i '_0.3i) 

0.66« 1.4) 


17.8!'^ 
~45 


I 42+<).40 
^•^^-0.38 

~6 


291/523 
2.3/3 


1.4 ±0.3 
«4.7 


Ch2 BB 


0.3« 1.1) 


~ 0.28^ 


R:i0.16 


2.4/3 


»0.90 


Ch2 mekal 


6+" 6 


~2.6 


0.7 ±0.2 


1.1/3 


RiO.54 



Note. — The uncertainties are given at the 68% confidence level for a 
single interesting parameter. For Ch2, the upper Hmits on nw value (at 68% 
confidence) are given in brackets, the lower limits are not constrained at the 
same confidence level. 

"NormaHzation for the PL model is the spectral flux at 1 keV in 10"* pho- 
tons cm"^ s"' keV"' . Normalization for the mekal model is the Emission 
Measure (EM) in 10^' cm"', scaled to a distance of 1 kpc. 

''Projected area of the emitting region for the BB model, in units of lO'' m", 
normalized to d = 1 kpc. 

■^BB temperature in keV. 

''We use the C statistic for Chi and the statistic for Ch2. 

■^Unabsorbed PL flux in the 0.5-8 keV band or bolometric BB flux, in units 
of 10"" ergs cm"^ s"' (corrected for vignetting and finite extraction aperture 
size). 



3. However, the short 8 ks exposure, the off-axis location, 
and the high EPIC background do not allow one to detect the 
faint extended component, while the compact PWN cannot be 
resolved from the pulsar because of the broad PSF of XMM- 
Newton. 

2.5. Optical-IR-mdio data 

To understand the nature of the X-ray sources Chi and Ch2, 
projected within the HESS source image and look for other 
sources that could be related to the TeV emission, we have 
examined the field at other wavelengths. We found no coun- 
terparts to Chi in the Two Micron All Sky Survey (2MASS; 
Skrutskie et al. 2006) or Digital Sky Survey (DSS2)^ catalogs, 
up to the limiting magnitudes Ks = 15.4, H = 16, J = 17.5, 
R = 19, and B = 21. The nearest optical/NIR source is a 
2MASS point source (J = 15.45 ±0.05, H = 13.10 ±0.05, 
K= 12.06 ±0.04) located at R.A. = 18^09"40.92^ decl. = 
-19°25'45.7". The 3.2" offset from the best-fit Chi position 
substantially exceeds the position uncertainty of « 0.4". Fig- 
ure 3 shows that Chi is projected very close to the extended 
Source A seen in the radio and IR images. 

The only NIR/optical source within the 5" radius of the Ch2 
position is a star, NOMAD 1 0705-0568334 in the NOMAD 
catalog (Zacharias et al. 2005), offset by only 0.9" from Ch2. 
Given the relatively large uncertainty of the Ch2 position (see 
§2.1), the positions of the star and Ch2 can be considered co- 
incident. The X-ray-to-optical flux ratio. Fx /Fy ^ 2 x lO"-', is 
typical for a K star with coronal X-ray emission (Maccacaro 
et al. 1988). The soft X-ray spectrum of Ch2 and the optical- 
NIR magnitudes of the stai- (B = 14.96, V = 14.66, R = 14.19, 
J= 13.35, H = 12.84, andK= 12.72) support such an interpre- 
tation. No Ch2 counterpart is seen in the NVSS 20 cm image 
or in the Spitzer 8 and 24 fim images. 

The NVSS 20 cm image (45" restoring beam size) and the 
25" resolution images by Brogan et al. (2004) show two bright 
compact radio sources (marked A and B in Fig. 3) projected 
near the HESS source center, with the 1 .4 GHz spectral fluxes 
of 0.3 Jy and 0.15 Jy for Source A and Source B, respec- 
tively. Within these sources, 6.7 GHz methanol masers have 
been detected (Pestalozzi et al. 2005): G10.95h-0.02 in Source 
A (peak flux 15 Jy at 6.7 GHz, Vlsr = 24± 1 km s"', near- 
distance 3.2 kpc) and Gl 1.03H-0.06 in Source B (peak flux 0.7 
Jy at 6.7 GHz, Vlsr = 20± 1 km s"', near-distance 2.9 kpc). 

Both A and B are also very bright far-IR sources, 
IRAS 18067-1927 and IRAS 18067-1921 in the IRAS Point 
Source Catalog v2.1 (IPAC Infrared Science Archive*"), with 
peak fluxes of 1292 and 513 Jy at 60 /im, respectively. The 
Spitzer IRAC 8 /im and MIPS 24 fim images (Fig. 3g,h; taken 
from the GLIMPSE survey^) show extended sources of irreg- 
ular shape at these positions, with characteristic sizes of ~ 1'. 
They are not seen in the DSS and 2MASS images, which sug- 
gests that they are intrinsically very cold (T ~ 70 K) and/or 
strongly absorbed molecular/dust complexes. 

The Spitzer images also reveal an extended source with 
an interesting morphology (shell-like in the 24 jim image), 
located 10' southwest of the HESS J1809 center (see 
Fig. 3g,li). The source has an extended radio counterpart 
(NVSS 180919-192904) clearly seen in the 20 cm image 
shown in the same figure. It has been detected by Helfand 
et al. (2006) in the Multi-Array Galactic Plane Survey (MAG- 

' see http://archive.eso.org/dss/dss 

* http://irsa.ipac.caltech.edu/ 

' http://www.astro.wisc.edu/sirtf/ 
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PIS) and proposed to be an SNR candidate GlO.8750+0.0875. 
It is possible that this extended source contributes to the non- 
uniform X-ray background on the SI chip (see §2.1 and Fig. 
3b). The SNR candidates Gll.03-0.05 and G11.18h-0.11 
(Brogan et al. 2004; see Fig. 3b) are too faint to be seen in 
the shallow NVSS images. 

3. DISCUSSION. 

The sub-arcsecond angular resolution of Chandra has al- 
lowed us to resolve the compact J 1 809 PWN, disentangle 
the pulsar and the extended emission components, and mea- 
sure their properties separately. Thanks to the very low ACIS 
background, we were also able to detect the faint large-scale 
emission surrounding the compact PWN. The observation has 
also provided serendipitous coverage of the central part of 
the extended TeV source HESS J1809 and revealed two X-ray 
sources bright enough to permit spectral measurements and 
accurate determination of their positions. In this section we 
discuss some implementations of our findings. 

3.1. TheJ1809PWN 
3.1.1. Luminosity and spectrum 

At the plausible distance of 3.5 kpc, the unabsorbed X-ray 
luminosity of the compact PWN, Lpwn ~ 4 x lO-'^ ergs s~', cor- 
responds to the X-ray efficiency, ripv/n=Lpvm/E ^ 2.2 x 10""^, 
similar to those of PWNe around younger and more powerful 
Vela-like pulsars (see Fig. 10). Some of the detected X-ray 
PWNe associated with older and less powerful pulsars have 
similar or higher X-ray efficiencies (e.g., rjpvm ~ 2 x 10"^ and 
4 x 10"* for J1509-5850 [r = 150 kyr] and B0355-h54 [r = 560 
kyr], respectively; Kargaltsev et al., in preparation), while 
others show much lower efficiencies (e.g., ?7pwn ~ 5 x 10"^ 
and 7 X lO"*" for PRS J1740H-1000 [r = 100 kyr] and Geminga 
[r = 340 kyr], respectively; Kargaltsev et al., in preparation; 
Pavlov et al. 2006). This indicates that jjpwn is not signifi- 
cantly correlated with r or E, at least for young and middle- 
aged pulsars (r < I Myr), perhaps because, in addition to E 
and r, it depends on other factors (e.g., the pulsar's speed and 
the angle between the spin and magnetic axes). 

The spectral slope of the J 1809 PWN, Fpwn = 1.4 ±0.1, is 
similar to those of the PWNe around Vela-like pulsars (listed 
in Table 2 of Kargaltsev et al. 2007a, hereafter KPG07a), ex- 
cept for two bow-shock PWNe with prominent tails, J1747- 
2958 (the Mouse; Gaensler et al. 2004) and B1757-24 (the 
Duck; Kaspi et al. 2001), which show softer spectra (F = 
2.0 ±0.2 and 2.5 ±0.3, respectively). Comparing the 11809 
spectrum with those of older PWNe, we see that it is similar 
to the spectra of relatively compact PWNe (e.g., F = 1 .4 ± 0.3 
for the bright part of the B0355H-54 PWN; McGowan et al. 
2006), but it is harder than the spectra of the extended tails 
in JI740H-I000 and J1509-5850 PWNe (F = 1.8 ±0.4 and 
2.2 ± 0.3, respectively; Kargaltsev et al., in preparation). It 
hints that PWN spectra are correlated with PWN morphology 
rather than with the pulsar age and spindown power. 

3.1.2. PWN morphology 

The unknown proper motion of J 1 809 and the faintness of 
the surrounding large-scale emission complicate the interpre- 
tation of the observed PWN, but its "cometary" appearance 
strongly suggests that the pulsar's motion plays a major role. 

3.1.2.1. Bow shock in a nearly isotropic pulsar wind? — 
The compact J1809 PWN is elongated approximately along 
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Fig . 10. — PWN luminosity versus pulsar spin-down power showing 1 1 809 
and 10 Vela-like pulsars observed with Chandra. The luminosities are esti- 
mated for the 0.5-8 keV band. The dash-dot lines are the lines of constant 
PWN efficiency, »jpwii. The error bars include the statistical uncertainties and 
the nominal 30% distance uncertainties, except for the Vela pulsar whose 
parallax has been measured (Dodson et al. 200i5a). 

the north-south direction, with the pulsar located much closer 
to its southern end (see Fig. I). Such a cometary morphology 
can be attributed to a bow shock created by the pulsar moving 
supersonically in the southern direction (P.A. » 194°). The 
apex of the termination shock (TS) of the pulsar wind is lo- 
cated at the distance 



Rh 



Efn 



-.1/2 



47rc(pamb + Pram) 



(1) 



ahead of the pulsar. At this distance, the pulsar wind pres- 
sure, Ph, = EfniAncr^)'^ (/o takes into account anisotropy 
of the pulsar wind), is balanced by the sum of the ambient 
pressure, pamb = pkT{iJ,mu)~^ = 1 .38 x 10"'^n^~' T4 ergs cm~^, 
and the ram pressure, /7ram = pv^ = 1-67 x 10~"'nv7 ergs cm"-' 
(T4 = T/10*K, V7 = v/10^cms~\ /i is the molecular weight, 
and n = p/niu is in units of cm"^). Assuming p^^^ ^ pamb 
(or ^ 1, where A4 = v/cs is the Mach number, c, = 

(5kT /Siimn)^^^ = 12/x"'/^r4 km s~^ is the sound speed in 

the ambient medium), we obtain Rh = l.7 x lO^'^n^^^^fQ^Vj^ 
cm. 

The shocked pulsar wind, observed in X-rays as a bow- 
shock PWN, is confined between the TS and the contact dis- 
continuity (CD) surface. For ^ I and a nearly isotropic 
preshock wind with a low magnetization parameter a, the TS 
acquires a bullet-like shape (Bucciantini et al. 2005, hereafter 
B05). The distance Rh between the pulsar and the bullet head 
is given by equation (I), while the CD surface head is at a dis- 
tance ~ 1.3/?/, from the pulsar. Assuming that the sharp rise 
in the brightness profile (Fig. 2) at about 3" from the pulsar 
(i.e. 1.6 X lO^^Js.s cm in the plane of the sky) corresponds to 
the CD head, we can estimate the pulsar velocity. 



140n"'/^J3 '5///^ sin / km s"' , 



(2) 



and the Mach number, M ~ l2n-^^% V^/Vn 4 
where i is the angle between the velocity vector and the line 
of sight. 

The observed projected length, ~ 12", and width, ~ 3", of 
the iimer PWN are close to the length ~ 7Rh ~ 15", and width. 
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^Rh'^ 2", of the TS bullet, predicted by the BOS models for 
(J a few X 10"^ The width of the outer PWN, ^ 20", can 
be interpreted as the diameter of the CD shell (approximately 
cylindrical behind the pulsar, with a diameter SRh ~ 18", 
according to BOS). 

The BOS models also predict a tail behind the back surface 
of the TS bullet, where the collimated shocked pulsar wind 
flows with subrelativistic velocities: 0.1c-0.3c in the inner 
channel with a cyhndrical radius ~ Rh 2" in our case), and 
up to 0.8c-0.9c in the outer channel, confined between the 
cylindrical surfaces with radii ~ Rh and ~ 4Rh. The lack of 
an extended tail in the ACIS images is puzzling; it could be 
due to a relatively low Mach number of the pulsar and insuf- 
ficient sensitivity of the short ACIS exposure. On the other 
hand, the shape of the faint thin feature northeast of inner 
PWN hints that it may be a bent extension of the inner PWN, 
which is difficult to explain if this is a strongly collimated, 
mildly relativistic flow just behind the back surface of the TS 
bullet (unless the flow is subject to kink instabilities). How- 
ever, statistical significance of this feature is too low to draw 
any definitive conclusions. 

3.1.2.2. Effects of pulsar wind anisotropy. — We should note 
that the BOS models assume an isotropic pulsar wind {fn =1). 
We know from observations of young pulsars, such as Crab 
and Vela, that the wind is not isotropic, but it is mostly con- 
fined to the equatorial plane. In addition, in young PWNe we 
often see jets along the pulsar's spin axis. These jets can be 
formed by polar outflows originating in the pulsar magneto- 
sphere or they can be created by tangential inflows just out- 
side the equatorial TS surface converging toward the spin axis 
(Komissarov & Lyubarsky 2004). We are unaware of PWN 
models that include both the pulsar wind anisotropy and the 
ram pressure effects. We can, however, expect that, in the case 
of an anisotropic wind, the PWN morphology, at least close 
to the pulsar, would depend on the orientation of the spin axis 
with respect to the pulsar's velocity. 

As Chandra observations have shown, the spin axis is ori- 
ented along the direction of pulsar's motion in a number of 
young pulsars (e.g., Ng & Romani 2004). In this case, we 
expect that the equatorial outflow would form a shell between 
the TS and the CD surface behind the pulsar, filled by a rel- 
ativistic plasma with a subrelativistic bulk flow velocity. The 
morphology of the X-ray emission from such a PWN would 
be generaUy similar to that in the case of isotropic wind, al- 
though we may expect to see a shell-like (rather than filled) 
PWN appearance at sufficiently high resolution. Such a struc- 
ture ("outer tails") is possibly seen in the Geminga PWN (Car- 
aveo et al. 2003; Pavlov et al. 2006). As for the jets coaligned 
with the pulsar velocity vector, the rear jet is expected to be 
seen along the axis of the shell filled by the shocked equato- 
rial wind behind the pulsar (at least if the jet originates from 
the pulsar magnetosphere), as observed in the Geminga PWN 
(Pavlov et al. 2006). The front jet is expected to pierce the 
head of the shell (as /n » 1 in the jet), and it can be seen 
ahead of the pulsar unless it is crushed by the ram pressure. 
The jets could be of quite different brightness and length be- 
cause of the Doppler boosting and different effective pres- 
sures in front and behind the pulsar 

In the case of J 1809, the high-resolution images shown in 
Figure 1 indicate the presence of a narrow structure (with a 
linear extent of ~ 3") just south of the pulsar, which might 
be interpreted as a front jet. Furthermore, the heavily binned 
low-resolution image (Fig. 3b) shows a faint structure extend- 



ing in the direction of the presumed pulsar's proper motion. 
Although tentative, the structure is consistent with the X-ray 
morphology seen in the ASCA GIS image (Fig. 3c). There- 
fore, one could speculate that these structures are connected 
to each other and could be parts of the front jet. In this in- 
terpretation, the rear jet may contribute to the inner PWN 
emission, and the above-mentioned bent extension of the in- 
ner PWN (§3.1.2.1) could represent the outer part of the rear 
jet. Although barely supported by the existing data, the jet 
hypothesis offers a way to explain the offset large-scale X-ray 
(and TeV) emission south of the pulsar, which cannot be as- 
sociated with the shocked pulsar wind under the assumption 
of wind isotropy. We speculate that this emission could be 
produced by particles supplied through the front jet (similar 
to the Vela PWN; Pavlov et al. 2003; Kargaltsev & Pavlov 
2004). However, to account for the luminous TeV emission 
(L~f ^2x lO'^d^^E), the front jet would have to carry a sub- 
stantial fraction of the pulsar's spindown power, unUke the 
northwest jet of the Vela PWN whose energy injection rate is 
only ~ 1Q~^E. 

We cannot also exclude a possibility that the spin axis of 
the J 1 809 pulsar, and hence the PWN jets, are substantially 
misaligned with the direction of pulsar motion. For instance, 
one could assume that the jets are nearly perpendicular to the 
velocity vector (i.e. this vector lies in the equatorial plane; a 
possible example is PSR B 1706^4 and its PWN; Romani et 
al. 200S). In this case, the jets would be bent backwards (or 
even destroyed) by the ram pressure of the oncoming ambi- 
ent medium, while the TS in the equatorial outflow (ring-like 
without ram pressure applied) would turn into an ellipse-like 
structure elongated in the direction of pulsar motion, with the 
pulsar displaced along the major axis toward the TS head. The 
observational appearance of the shocked pulsar wind would 
then depend on the inchnation of the orbital plane to the line 
of sight, but at most inclinations we would see an elongated 
structure with the pulsar shifted from the center of the PWN in 
the direction of the velocity vector, in a qualitative agreement 
with the observed shape of the inner PWN. Of course, such an 
interpretation would imply a higher pulsar velocity than that 
estimated above for a nearly isotropic pulsar wind (see eq. 
[2]). On the other hand, in such a geometry the wind could 
propagate to larger distances ahead of the pulsar because the 
CD could be easier destroyed by various instabilities, and 
the pulsar wind would be mixed with the shocked ambient 
medium. We note, however, that this picture can hardly ex- 
plain the large-scale X-ray emission south-southwest of the 
pulsar, which, in this case, could be attributed to the host SNR 
or the crushed relic PWN (see §3.3 and §3.4.2, respectively). 

Overall, we can conclude that the observed PWN morphol- 
ogy is generally consistent with the assumption of the su- 
personic pulsar motion, which implies a low ambient pres- 
sure (e.g., /5an,b ^ 2 X 10~'*'/^^'^3 5 sin^ ' ergs cm"-' for a nearly 
isotropic pulsar wind). On the other hand, modeling of mag- 
netized anisotropic winds from fast-moving pulsars, deeper 
X-ray observations, and proper motion measurements are 
needed to firmly establish the nature of the J1809 PWN and 
infer its properties quantitatively. 

3.2. The pulsar 

The spectrum of the J1809 pulsar is less certain than that of 
the PWN. Although the one-component PL model formally 
fits the spectrum, the fit yields a rather large F ~ 2.6-3.2 
and suggests a smaUer nu than the one obtained from the 
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PL fit to a better quality PWN spectrum. An alternative 
description of tlie pulsar spectrum is provided by the two- 
component PL+BB model that is often used to fit the spectra 
of young and middle-aged pulsars (see, e.g., KPG07a and ref- 
erences therein). The PL component of the PL-hBB fit gives 
Fpsr ~ 1.2 ±0.6, similar to Fp^n and to the spectral slopes 
of other pulsars of similar ages (see e.g., KPG07a). In any 

•-1/2 

case, the slope is much softer than Fpsr = 2.1 — 2.9£3g « 0.1 
predicted by the Gotthelf 's (2003) correlation. The X-ray effi- 
ciency of the pulsar, ?7psr = L^si/E ^ 2 x lO^^tif j in the 0.5-8 
keV band, is also not unusual for Vela-hke pulsars, as well as 
the ratio Lpwn/^rsr ~ 10 (see KPG07a). 

The possible thermal component of the PLh-BB fit is poorly 
constrained, not only because of the scarce statistics but also 
because the soft thermal radiation is strongly absorbed by the 
ISM. The BB temperature, T ~ 1.7-2.3 MK, emitting area 
A ~ lO^-lO'' m^, and bolometric luminosity, L^^ ~ (0.6- 

1.6) X 10^^ ergs s"', are similar to those found from the 
PLh-BB fits for Vela-like pulsars. The temperature is a fac- 
tor of two higher than the NS surface temperatures predicted 
by standard NS cooling models for the 50 kyr age (e.g., 
Yakovlev & Pethick 2004), and the corresponding emitting 
area is smaller than the NS surface area. However, the actual 
spectrum of the NS thermal radiation can differ substantially 
from the BB model. In particular, fitting the spectra with hy- 
drogen atmosphere models (Pavlov et al. 1995) yields lower 
effective temperatures and larger emitting areas (see Pavlov et 
al. 2001a for the specific example of the Vela pulsar). Unfor- 
tunately, the quality of the data does not warrant fits with more 
complicated atmosphere models. Unlike the temperature and 
the area, the bolometric luminosity is not so sensitive to the 
presence and properties of the NS atmosphere. The compari- 
son with the NS cooling models shows that the J1809's bolo- 
metric luminosity is a factor of ^ 10 lower than predicted by 
the "basic" theoretical cooling curve (Mns = 1.3M0, no su- 
perfluidity); it is consistent with the cooUng curves for heavier 
NSs (e.g., Mns = 1-5-1.6M0) for various superfluidity mod- 
els (Yakovlev & Pethick 2004). 

3.3. HostSNR 

With the velocity given by Equation (2), the pulsar would 

have traveled a distance of ~ ln~^l^dl\fll^ sin i pc during the 
time equal to its spindown age, r = 5 i kyr. This corresponds 

to a displacement of ~ l'rr^l^dl\fl[^svc^i in the plane of 

the sky and, for n"'/^c/3 5/j/^sin^/ ~ 0.6, places the sky pro- 
jection of the pulsar's birthplace close to the apparent center 
of die radio SNR Gl l.lS-nO.l 1 (see Fig. 3b), suggesting that 
J1809 and G11.18h-0.11 were created by the same supernova 
explosion. This conjecture implies that the pulsar has over- 
taken the SNR shell and left the high-pressure SNR interiors, 
which is consistent with the assumption /?a,nb <^ Pnm, used to 
obtain the pulsar velocity estimate. This could explain why 
J 1809 looks so different from those Vela-hke PWNe fliat do 
not show cometary morphology because they are moving sub- 
sonically in high-pressure SNR interiors (see KPG07a for ex- 
amples). However, if G11.18h-0.11 is at li « 3.5 kpc (the pre- 
sumed distance to J 1809), then its size, D « ^d^.s pc, would 
be surprisingly small for a 50 kyr old SNR, unless it is ex- 
panding in an unusually dense environment. To get more con- 
sistent SNR age and size, we have to assume a larger distance 
and a smaller pulsar's true age, which would require a lower 
ambient density around J1809 to match the pulsar's birth- 



place with the center of G 11.18+0.11, resulting in a higher 
inferred pulsar velocity and Mach number. For instance, if 
we assume d = 5 kpc (i.e., Z) ~ 1 1 pc) and the true age = 20 
kyr, then n ~ O.l/o sin'' i cm~^ is required, which corresponds 
to V - 300(sin/r' km s"', - 25T^^^^ fi^^^isiniTK Thus, 
given the uncertainty of the pulsar's true age and distance, we 
cannot rule out the possibility that Gil. 18+0.01 is the host 
SNRforJ1809. 

It is also possible that G 11.18+0.01 is just a background 
SNR accidentally projected near the pulsar. In this case, the 
real host SNR might be associated with the extended emis- 
sion around J1809 seen in the Chandra and ASCA images 
(Fig. 1, bottom right and Fig. 3c). The faintness of the large- 
scale X-ray emission does not allow one to determine its ori- 
gin (thermal or nonthermal) and, therefore, thermal plasma 
emission from the SNR interior cannot be excluded until bet- 
ter quahty data are obtained (see §2.2.1). However, the size of 
this putative SNR is too smaU, and, more importantly, the ob- 
served PWN morphology does not look consistent with this 
hypothesis. Therefore, it seems more plausible that the ex- 
tended emission is synchrotron radiation related to the PWN 
(see §3.1.2.2) rather than thermal emission from hot gas in 
SNR interior. Measuring the X-ray spectrum of the extended 
emission in a deeper observation would distinguish between 
these possibilities. 

One could also speculate that the host SNR is not seen be- 
cause its size is larger than the field-of-view of our observa- 
tion (> 20 pc at c/ = 3.5 kpc), and the SNR interior has be- 
come cold enough (r < 10^ K) to provide M'> \ and be 
undetectable in X-rays. This option remains quite viable if 
the pulsar's true age is close to (or exceeds) its spindown age. 

3.4. Origin of HESS Jl 809 

TeV radiation can be produced by the inverse Compton 
scattering (ICS) of low-frequency radiation (e.g., the cos- 
mic microwave background radiation [CMBR]) off relativistic 
electrons. Alternatively, it can be generated by the tt" — > 7+7 
decay, the tt" mesons being produced when relativistic pro- 
tons interact with the ambient matter. Therefore, to under- 
stand the nature of a TeV source, one should identify the 
source of relativistic particles and the target with which these 
particles interact. 

3.4.1. Possible targets where the TeV radiation is produced 

The omnipresent target for relativistic electrons producing 
TeV photons by the ICS is the CMBR, with the energy den- 
sity t/cMBR = 0.26 eV cm~^. To produce TeV photons with 
energy E-y by upscattering the CMBR photons with energy 
e ^ 3kT ^ 4 X 10~^ eV, electrons with the Lorentz factor 
7 ~ 5 X lO^EjJ^ [i.e. Ee ~ 25EjJ^ TeV] are required, where 
£'TeV= £'7/(1 TeV). As 76 ~ 20-100 keV-Cm^c^ for ~ 1- 
20 TeV, the ICS occurs in the Thompson regime. The same 
relativistic electrons produce synchrotron photons with en- 
ergy 

£sy„ - I'h^cyc ~ 57?B-5 eV ~ QA-j^^f^^^kcY , (3) 

where 77 = 7/10^ andfi-5 = B/ilOpG). 

Close to the Galactic plane, where HESS Jl 809 is situ- 
ated, a factor of a few higher radiation energy density can 
be provided by IR emission from interstellar dust and Galac- 
tic starlight. For instance, the models of intersteUar radiation 
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field by Strong et al. (2000) give the energy densities of 0.6 
and 2.7 eV cm"^^ for these two components, respectively, at 
the galactocentric distance of ~ 4 kpc. Since the ICS pro- 
ceeds in the Klein-Nishina (K.-N.) regime for e > l^xev 
(i.e., the K.-N. effects become important at e > a few x 10~^ 
eV for the high-energy end of the TeV photon spectrum), we 
should use a more general formula for estimating the electron 
Lorentz factor: 



7~ 10^ 



fixeV + (£^TeV + ^Tev/ eev) 



1/2- 



(4) 



where Cev = e/(leV). The electrons with such a 
Lorentz factor generate synchrotron photons with energy 



■^syn 



0.05 



,1/2 



B-5 eV. For in- 



EjeV + (^TeV + ^Tev/ Cev) 

stance, electrons with 7 ^ 2 x 10^ (i.e. Eg ~ 40 TeV) pro- 
duce photons with = 10 TeV by upscattering background 
starlight photons with a typical energy e ~ 1 eV, and the same 
electrons produce synchrotron radiation in a far-UV range, 
Esyn ~ 6 eV, in an interstellar field of 3 /xG. Although the radi- 
ation energy density of the starlight can substantially exceed 
that of the CMBR, the emissivity is reduced by the smaller 
Compton cross section in the K.-N. regime and smaller pho- 
ton number density. 

The radiation energy density and photon number density 
(hence TeV emissivity) can be enhanced in the vicinity of very 
bright sources. The radiation energy density at an angular 
distance 6 from a source of radiation can be estimated asU = 
F /{cO'^), where F is the source energy flux observed at Earth; 
it exceed the ambient radiation energy density, t/amb> at 



> (F/cU,^b) 



1/2 



(5) 



The multiwavelength data show two very bright radio-lR 
objects near the center of HESS J 1809 - Source A and Source 
B (see Fig. 3), which are likely molecular-dust complexes in 
star-forming regions. An approximate integration of the spec- 
tral flux of Source A gives its total flux Fa ^ 10"^ ergs cm"^ 
s~' (i.e. the luminosity La ^ 1.5 x lO^^d^^ ergs s"'). This 
means that the radiation energy density near this source ex- 
ceeds Usanb at6 < 5'(f/amb/l eVcm"^)"'/-^. For instance, if the 
TeV emission is produced by the ICS of IR photons with ener- 
gies around e ^ 0.02 eV (which corresponds to the maximum 
spectral flux of Source A), and U-^iab ^ OA eV cm"^ in this 
energy range, then U > U-^mh at < 15'. This value is close 
to the observed radius of HESS J 1809, which suggests that 
the sphere of ~ 10 pc of the enhanced radiation field around 
Source A might be the target where the TeV emission is pro- 
duced by the ICS off the relativistic electrons supplied by a 
nearby source [Lorentz factors 7 ^ (0.5-6) x 10^ are required 
to upscatter the IR photons to Ej = 0.5-20 TeV]. Similar es- 
timates hold for Source B, whose flux is a factor of 2.5 lower, 
and the size 6* is a factor of 1 .6 smaller, than those for Source 
A. The hypothesis that the "photon sphere" of Source A or 
Source B is the site where the TeV emission is produced via 
ICS could be verified by detection of the accompanying syn- 
chrotron radiation. However, most of the energy range of the 
corresponding synchrotron emission, Zisyn ^ (1-200)B_5 eV, 
is subject to the strong interstellar absorption. 

In principle, the TeV emission could be initiated by high- 
energy nucleons suppUed by some source (e.g., a pulsar), 
which produce decaying tt" mesons in collisions with nucle- 
ons of the ambient medium (e.g.. Horns et al. 2006). The 
target in this case would the nucleonic component of the 
circumpulsar medium, and we should expect enhanced TeV 



emission in dense, cold clouds, if there are such clouds close 
to the source of high-energy nucleons. The lack of an ex- 
tended IR or radio counterpart of a size similar to that of 
HESS Jl 809 (which would indicate the presence of a large 
molecular cloud) suggests that the TeV source is not associ- 
ated with a localized target for high-energy nucleons. 

3.4.2. Possible sources of relativistic particles 

In the case of HESS J1809, the natural source of relativistic 
electrons (and possibly protons) is the J 1809 pulsar. In ad- 
dition to that, there are the Chi and Ch2 X-ray sources and 
the radio SNR Gl 1.03-0.05, all projected close the center of 
HESS J809 (see Fig. 3). Finally, relativistic particles might be 
produced by acceleration mechanisms in the forward shock of 
the putative host SNR of the J1809 pulsar. 

The soft X-ray spectrum of Ch2 and the positional coinci- 
dence with a field star suggest that the X-ray emission comes 
from an active stellar corona. Therefore, we conclude that 
Ch2 is not related to HESS J1809. 

If the sheU-like radio SNR Gl 1.03-0.05 were a power- 
ful source of relativistic electrons, its synchrotron radiation 
would have been seen in the X-ray range. Since no X-rays 
at the SNR location are detected in our Chandra observation, 
this SNR is not a viable candidate for the source of relativistic 
particles fliat powers HESS J 1809. 

As we have discussed in §3.1.2 and §3.3, the faint, large- 
scale X-ray emission south of the J 1809 pulsar is likely not 
the host SNR of this pulsar, but it is rather related to the J 1809 
PWN. Moreover, TeV 7-ray emission from an SNR shock is 
usually associated with the SNR shell, not the interior. Cur- 
rent radio. X-ray and 7-ray data provide no evidence of such 
a shell around J 1809. Therefore, even if the large-scale X-ray 
emission belongs to the host SNR, HESS J 1809 is unlikely to 
be powered by particles accelerated in the forward shock of 
this SNR. 

Among the X-ray sources detected with Chandra, Chi is 
the closest to the projected HESS J1809 center (Fig. 3). The 
lack of an IR-optical counterpart and the large lower limit 
on the X-ray-to-optical flux ratio, Fx/Fopt ^ 3, mean that the 
X-ray emission from Chi comes not from a usual field star. 
On the other hand, the absorbing hydrogen column density, 
nH.22 = 1.2 ±0.4, is close to that of the J1809 PWN, which 
suggests that Chi is a Galactic object, not an AGN observed 
through the Galactic plane. The slope of the Chi spectrum, 
r = 1 .4 ± 0.4, suggests that it may be either an X-ray binary 
or a remote PWN unresolved because of the far off-axis lo- 
cation. With the optical extinction estimated from the mea- 
sured tin value (e.g.. Ay ~ 4-9), the limits on the optical-NlR 
magnitudes (see §2.5) virtually exclude a high mass X-ray 
binary. Although a low mass X-ray binary (LMXB) in a qui- 
escent state could have an X-ray luminosity comparable to 
that of Chi, Lchi 3 x 10^2(^/4 ^^^^^2 g^g^ ^-i^ LMXBs are 
not known to be sources of TeV emission. If Chi is a young 
pulsar with a PWN, then it might provide relativistic particles 
needed for generating the TeV emission. (We should note, 
however, that the TeV-to-X-ray flux ratio, ~ 150, would be 
much higher than the values, ^ 0.01-3, inferred for most of 
the other TeV PWNe, except for the B 1800-21, for which 
the ratio is about 100; see Table 2 in Kargaltsev et al. 2007b, 
hereafter KPG07b, and references therein). To assess the like- 
lihood of Chi being the source of relativistic particles for 
HESSJ1809, a deeper Chandra observation is needed, with 
Chi imaged close to the telescope's optical axis. 

Another plausible physical counterpart to HESS Jl 809 is 
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the J 1809 pulsar/PWN. To date, young pulsars have been 
found in the vicinity of ^ 10 extended TeV sources (e.g., de 
Jager 2006; Gallant 2007), and the hkehhood of this happen- 
ing by chance is very low (see, e.g., KPG07b). The ratio 
of the TeV flux of HESS J1809 to the X-ray flux die J1809 
PWN (including the outer component of the compact PWN) 
is about 100, higher than those observed in most TeV ple- 
rions (KPG07b). However, if one adds in the X-ray flux of 
the large-scale extended emission (§2.2.1), then the ratio de- 
creases down to ~ 20 (or even lower if one includes the diffuse 
emission of a larger extent, seen by ASCA; see Fig. 3c). 

The center of the brightest part of HESS J 1809 is offset 
from J 1809 by « 8', almost in the direction of the presumed 
pulsar motion. Such an offset is comparable to the 10'-20' 
offsets found in most of TeV plerions (including the most 
secure associations PSRB0833-45/HESS J0835-455 [VelaJ 
and PSRB1823-13/HESS J1825-137). An obvious explana- 
tion for such an offset is that the photons for the ICS are pro- 
vided by a compact source (e.g.. Source A) whose center lies 
at some distance from the pulsar (> Sds s pc for Source A). 
Also, the offset might be due to the anisotropic supply of rel- 
ativistic particles through the putative front jet (see §3.1.2). 

If the TeV emission is produced by the ICS of an approxi- 
mately uniform background radiation, such as the CMBR or 
Galactic starlight, then the origin of the offsets and the asym- 
metries of the extended X-ray and TeV PWN components 
could be attributed to the reverse SNR shock that had propa- 
gated through the nonhomogeneous SNR interior and reached 
one side of the PWN sooner than the other side, crushing the 
PWN and sweeping its contents along (Blondin et al. 2001). 
One can assume that the TeV source HESS J1809 is powered 
by this relic crushed PWN, whose relativistic electrons have 
not lost all their energy to synchrotron radiation because of, 
e.g., a lower magnetic field, while the compact PWN in the 
pulsar vicinity is created by fresh electrons, recently injected 
from the pulsar magnetosphere. As the electrons responsi- 
ble for the TeV emission have been accumulating during a 
substantial fraction of the pulsar's lifetime, and the pulsar's 
spindown power was higher in the past, this hypothesis could 
explain the uncomfortably large ratio, ~ 0.01^1 5, of the TeV 
luminosity to the current spindown power. (A similar expla- 
nation has been suggested by Aharonian et al. 2006 for the 
PSRBl 823-1 3/HESS J1825-137 association.) This hypoth- 
esis can be tested with deep radio observations, which should 
be able to detect the synchrotron emission from the rehc elec- 
trons of the crushed PWN. 

In the above picture, the asymmetric extended X-ray emis- 
sion south of the J 1809 pulsar could also be associated with 
the crushed PWN and considered as an X-ray (synchrotron) 
counterpart of the TeV (IC) source HESS J1809. The smaU 
size of the extended X-ray emission, in comparison with the 
size of the TeV source, and its offset from the center of 
HESS J 1809 could be attributed to the lower sensitivity of the 
ACIS chip S2 (where the central part of HESS J1809 is im- 
aged - see Fig. 3a) or to a lower magnetic field in the central 
part. A deeper X-ray observation could distinguish between 
these two possibilities. 

The large-scale morphology of the X-ray emission around 
J1809 is remarkably similar to that around PSR B1823-13, 
which powers the X-ray PWN G18.0-0.7 (Lx ~ 3 x lO^^ ergs 
s~\ angular size > 5'; Gaensler et al. 2003). In addition to 
the extended low-surface-brightness component, the B1823- 
13 PWN has a much more compact (5"-10") brighter core, 
resolved by Chandra (Teter et al. 2004). The TeV emission 



from HESS 1825-137 covers an area much larger than the 
X-ray PWN area, extending up to 1 ° southward from the pul- 
sar (Aharonian et al. 2006). However, just as in the case of 
PSR J1809/HESS J1809, both the TeV and the low-surface- 
brightness X-ray emission have similar shapes, and they are 
offset in the same direction with respect to the pulsar position. 

Finally, we would like to point out that there are two other 
young pulsars, J181 1-1925 and J1809-1943, at distances of 
« 24' and « 19', respectively, from the HESS J 1809 cen- 
ter. PSR Jl 8 11-1925 can not be related to HESS J 1809 since 
it has been associated with Gl 1.2-0.3 (Kaspi et al. 2001), 
whose diameter, w 4', is substantially smaller than the offset 
from HESS J1809. PSR J1809-1943 is the radio counterpart 
of the transient Anomalous X-ray Pulsar (AXP) XTE J1810- 
197 (Halpem et al. 2005, and references therein). Its associ- 
ation with HESS J1809 does not look very plausible because 
of the large spatial offset and because its properties are quite 
different from those of the young pulsars found in the vicinity 
of other TeV sources (no TeV emission from AXPs have been 
reported so far). 

Thus, based on the current data, only two objects, the J1809 
pulsar/PWN and, less likely, Chi, are plausible candidates for 
the source of relativistic particles powering HESS J 1809. In 
both cases, the TeV emission is generated by ICS of either IR 
photons from Source A (or Source B) or CMBR photons. To 
discriminate between different possibilities, the nature of Chl 
should be established in deep, high-resolution Chandra ob- 
servations, and the extended X-ray PWN emission should be 
studied in a deep XMM-Newton exposure. The synchrotron 
emission from the alleged relic PWN could be detected in 
deep radio observations. 

4. CONCLUSION 

We have detected the X-ray emission from PSR B 1809-19 
and its synchrotron nebula. The X-ray efficiency and spec- 
trum of the PWN are similar to those of many other compact 
PWNe, both younger and older than J1809. The cometary 
shape of the compact PWN suggests that the pulsar is moving 
supersonically, but no extended tail is detected behind the pul- 
sar, perhaps because the Mach number is not large enough or 
the exposure time is too short. Our analysis of the PWN mor- 
phology suggests an anisotropic pulsar wind outflow, possibly 
including jets oriented along the direction of pulsar motion. 
To establish the nature of the compact PWN unambiguously, 
proper motion measurements and a deeper Chandra observa- 
tion are required. 

The compact PWN is immersed in an extended emission 
of lower surface brightness. This extended emission is offset 
with respect to the pulsar in the direction opposite to that of 
the compact PWN (i.e. in the direction of the alleged proper 
motion). If the extended emission is indeed powered by 
the supersonically moving pulsar, then the radiating particles 
could be supplied through the front jet of the compact PWN 
or brought by the northern part of the reverse SNR shock that 
overtook the pulsar moving southward. 

The spectrum of the pulsar can be described by a two- 
component BBh-PL model. For reasonable hh values, the pa- 
rameters of this model and the corresponding component lu- 
minosities resemble those of Vela-like pulsars. 

The J 1809 pulsar and its PWN are located in the vicin- 
ity of tiie unidentified extended TeV source HESS J1809- 
193, which can be powered by ultrarelativistic electrons 
accelerated in the J 1809 PWN. In addition to J 1809, we 
found another X-ray source, Chl, within the central part of 
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HESS J 1809, which might also be an unresolved pulsar/PWN. 
Although the true nature of Chi remains elusive, we cannot 
rule out the possibility that Chi is an alternative X-ray coun- 
terpart of HESS J1809. Whatever of the two objects, J1809 
or Chi, supply the relativistic electrons, the TeV emission is 
likely produced by the ICS of either CMBR or IR photons 
from bright IRAS sources (Source A or Source B) projected 
near HESS J1809. 
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